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ABSTRACT 


Equilibrium  phase  compositions  and  refractive  indices 
were  determined  for  the  carbonyl  sulfide-propane  system  using  a 
windowed  cell  of  variable  volume.  Measurements  were  made  at 
temperatures  of  20°,  60°,  100°,  140°,  and  180°F  and  at  a  series  of 
pressures  ranging  from  the  vapor  pressure  of  propane  to  the  vapor 
pressure  of  carbonyl  sulfide  at  each  temperature. 

Pressure-composition  diagrams  for  each  isotherm  showed  an 
S-shaped  configuration  for  both  the  saturated  vapor  and  saturated 
liquid  loci.  Contrary  to  expectation,  no  azeotropes  were  encountered, 
although  a  tendency  toward  azeotropic  behavior  existed  at  the  high 
concentration  end  for  carbonyl  sulfide. 

The  molar  volumes  of  the  co-existing  phases  were  calculated 
from  the  measured  refractive  indices  and  equilibrium  phase 
composition  data  using  the  Lorentz-Lorenz  relationship.  The  pure 
component  refractivi ti es  were  assumed  to  be  additive  on  a  molar  basis. 

Early  in  the  study,  it  became  apparent  that  significant 
errors  existed  in  the  vapor  pressures  and  critical  properties  reported 
for  carbonyl  sulfide.  Accordingly,  the  vapor  pressure  and  the 
refractive  index  of  both  liquid  and  vapor  for  pure  carbonyl  sulfide 
were  measured.  The  new  vapor  pressure  data  are  reported  from  20°F 
to  the  critical  temperature.  The  critical  pressure  was  estimated  to 
be  960  psia  and  the  critical  temperature  to  be  235° F.  This  compares 
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to  earlier  calculated  values  of  897  psia  and  22]°F  respectively. 

The  experimental  vapor  pressure  data  for  carbonyl  sulfide 
were  used  to  determine  a  value  of  0.056  for  Pitzer's  accentric 
factor  to.  This  value,  together  with  the  critical  properties,  made 
it  possible  to  estimate  a  value  of  0.0044  for  the  k .  .  interaction 

'  J 

parameter  for  mixtures  of  propane  and  carbonyl  sulfide. 

The  equilibrium  data  obtained  at  100°F  were  used  to 
estimate  the  constants  needed  in  the  Chueh-Prausni tz  computer 
program  developed  for  predicting  the  phase  behavior  of  multi- 
component  systems.  These  constants  were  used  to  predict  the  data 
of  this  investigation  at  the  remaining  four  isotherms.  Close 
agreement  between  the  experimental  and  predicted  equilibrium  ratios 
resulted.  A  thermodynamic  consistency  test  was  made  on  the  data 
using  the  technique  presented  by  Chueh,  Muirbrook,and  Prausnitz. 

The  inconsistency  ranged  from  5.7  to  22.7  percent  as  the  temperature 
changed  from  20°  to  180°F. 
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I  -  INTRODUCTION 

Vapor-liquid  equilibrium  relationships  are  essential  for 
the  solution  of  many  process  engineering  problems.  In  the  design 
of  distillation  towers,  absorbers,  strippers,  and  other  contacting 
equipment,  a  knowledge  of  the  composition  of  the  stream  being 
treated  and  of  the  relative  amounts  of  the  co-existing  phases  under 
equilibrium  conditions  is  needed. 

The  data  required  for  these  relationships  are  often 
determined  experimental ly.  Over  the  years,  a  considerable  amount 
of  vapor-liquid  equilibrium  data  have  been  accumulated  by  using  a 
variety  of  techniques.  Nevertheless,  there  are  still  a  large 
number  of  industrially  important  systems  whose  equilibrium  relation¬ 
ships  have  not  yet  been  studied.  Among  these  are  the  carbonyl 
sulfide-hydrocarbon  systems. 

Carbonyl  sulfide  is  reported  to  be  a  common  contaminant 

of  refinery  and  synthesis  gases^.  It  is  believed  to  form  during 

thermal  or  catalytic  cracking  and  reforming  operations  and  it  exists 

in  concentrations  which  vary  from  10  to  20  ppm.  Like  other  sulfur 

contaminants  in  hydrocarbon  fluids,  its  presence  is  detrimental  to 

the  processing  of  the  fluids  or  marketing  of  products  derived  from 
2 

them  .  Its  vapor  pressure  characteristics  and  critical  properties 
are  very  close  to  that  of  propane,  thereby  posing  problems  in  the 

3 

purification  processes  of  the  lower-boiling  products  .  Previous 
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studies  on  its  separation  showed  that  it  is  not  satisfactorily 
removed  by  the  conventional  amine  scrubbing  processes  employed  by 
the  industry  for  the  removal  of  H^S,  mercaptans,  and  other  sulfur 
bearing  compounds.  Instead,  it  reacts  with  amine  solutions  to  form 
stable  compounds  by  a  complex  re-arrangement  so  that  the  de¬ 
activated  amine  solutions  can  not  be  recovered  by  regeneration  or 

1  2 

revivification  processes  *  . 

An  experimental  program  to  study  the  phase  behaviour  of 
the  carbonyl  sulfide-propane  system  was  thus  developed.  The  data 
obtained  may  then  be  used  to  estimate  the  binary  interaction  co¬ 
efficient  k..  and  other  constants  required  for  the  analytical 

modelling  of  multicomponent  phase  behaviour  according  to  the  Chueh- 

4  5 

Prausmtz  and  the  Soave-Redl ich-Kwong  methods. 
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II  «  THEORY 


A)  Refraction 

When  light  travels  from  one  medium  to  another,  it 
undergoes  a  change  of  velocity,  and  if  the  angle  of  incidence  is 
not  90°,  there  will  be  a  change  in  direction.  This  phenomenon 
is  commonly  known  as  refraction. 

The  law  of  refraction,  widely  referred  to  as  Snell's 
law,  has  been  attributed  to  W.  Snell  though  other  scholars 
claimed  that  it  was  Descartes  who  made  the  first  discovery  in  1637. 

In  any  event,  this  law  has  been  the  basis  for  the  determination  of 
refractive  indices. 

Figure  1  demonstrates  a  ray  diagram  from  which  an 
expression  for  the  refractive  index  may  be  deduced.  From  Snell's 
law,  it  follows  that: 

n-j  siny  =  sing 
sing  =  n^  sin6 

n^  =  n-j  siny/sin<5  (1 ) 

where  n-j  is  the  refractive  index  of  the  ambient  air,  n£  the  refractive 
index  of  the  glass  window,  and  n^  the  refractive  index  of  the  fluid 
being  measured.  The  angle  6  is  the  prism  angle  and  y  is  the  angle 
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FIG.  1  SCHEMATIC  REPRESENTATION  OF  THE  ANGLE  MEASUREMENT  INVOLVED 
IN  THE  DETERMINATION  OF  THE  REFRACTIVE  INDEX 


.  .  . 


. 


' 


5 


measured  by  autocol 1 imation .  For  a  light  source  of  given  wavelength, 
n-j  is  known,  6  can  be  measured  if  n^n-j,  that  is,  if  the  cell 
contains  ambient  air,  and  hence,  n^  can  easily  be  determined. 


B)  Density  Measurement  From  Refractive  Index 

Newton'7  in  1704  introduced  the  concept  of  refraction  in 
terms  of  refractive  index  and  density.  His  formulation  shows  that: 


n-1 

d 


r^  =  constant 


(2) 


O 

About  a  century  later,  Laplace  in  1806  postulated  that-  the 
proportionality  constant  r-j  is  invariant  with  pressure  and  temp¬ 
erature.  Later  investigations^’ ^  however  ,  indicated  that 
Laplace's  conclusion  is  not  always  valid  but  variation  of  r-j  is  so 
small  that  it  may  be  assumed  to  be  constant  without  introducing 
serious  error. 


Since  then,  numerous  formulae  relating  density  and  re- 

12 

tractive  index  have  been  proposed.  Besserer  has  reported  a 
comprehensive  review  of  the  more  relevant  refractivity  formulae. 
Among  these  correlations,  the  one  proposed  by  Lorentz^  and  Lorenz^ 
is  widely  accepted  in  optical  physics  mainly  because  it  can  be 
related  to  the  electronic  polarizability.  The  relationship  is: 


n2-l 

n2+2 


-  =  R 


LL 


M 

P 


(3) 
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where  R^  is  the  molar  refractivity .  Batsanov^  outlined  the 
derivation  of  this  formula  from  electromagnetic  theory. 

The  relationship  of  refractive  index  and  density  was 

originally  tested  on  pure  substances.  Subsequently,  scientists 

established  the  effect  of  composition  on  the  refractivity  of 

1  17 

mixtures.  The  work  of  Biot  and  Arago  ’  ;  Smith,  Engel,  and 

18  19  20  21 

Wilson  ;  Young  and  Finn  ;  Keilich  ;  and  Bloom  and  Rhodes  showed 

that  pure  component  ref ractivities  are  nearly  additive  in  their 

molar  ratios  to  give  refract! vi ties  of  mixtures.  The  deviations 

from  this  assumption  are  small  (about  1  to  2%)  and  within  the  accuracy 

of  most  experimental  techniques.  Thus,  the  assumption  that  mixture 

refractivities  can  be  expressed  by  a  molar  sum  of  the  pure  component 

refractivities  is  justifiable. 

C)  Vapor-Liquid  Equilibrium 

Generally  speaking,  a  system  is  said  to  exist  in 
equilibrium  if  no  change  in  the  condition  of  the  system  occurs  with 
time.  These  changes  are  brought  about  by  resistances  and  by  driving 
forces.  Some  of  these  forces  are  attributed  to  pressure  and/or  temp¬ 
erature  gradients  and  chemical  potentials. 

In  reality  however,  a  true  state  of  equilibrium  is  probably 
never  reached.  There  is  a  continual  variation  in  our  surroundings 
caused  by  resistances  which  have  a  tendency  to  diminish  with  time. 


is  *  '  >  ■  -  • 
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Nevertheless,  as  equilibrium  is  approached  the  rates  of  change 

become  very  small  so  as  to  become  immeasurable  by  available  means. 

For  practical  purposes  this  state  might  be  considered  as  a  state  of 

equilibrium.  In  many  engineering  problems  the  assumption  of 

equilibrium  is  justified  when  the  results  assuming  a  state  of 

22 

equilibrium  are  of  satisfactory  accuracy  . 

The  condition  of  equilibrium  can  be  expressed  in  terms 
of  thermodynami c  variables.  Combining  the  first  and  second  laws 
of  thermodynamics,  it  can  be  shown  that: 

<d  “system  +  Pd  System  ’  Td  Ssystem>  =  0  <4> 

Various  restraints  may  be  put  on  the  above  expression  to  produce 
alternative  conditions;  for  example: 

(d  “system)  „  =  0  <5> 

system,  vsystem 

A1 ternati vely,  other  pairs  of  properties  may  be  held  constant.  The 
most  useful  result  comes  from  fixing  T  and  P.  In  so  doing, 

d(U  .  +  P  V  .  -TS  +  )  =0  (6) 

v  system  system  system^  p  ' 

The  term  inside  the  bracket  is  commonly  known  as  the  Gibbs  Free 
Energy  (G).  Therefore,  for  a  system  to  be  in  equilibrium: 


■ 


d^system^  p  0 


(7) 


The  application  of  this  criterion  to  a  multiphase  closed  system  gives 
the  generalised  condition  of  phase  equilibrium: 

i  ii  iii 

ui  =  ui  =  vi  ....  (8) 


where  u. 1 s  are  the  chemical  potentials  of  i  in  each  phase. 

A1  ternatively,  equation  (8)  can  be  transformed  to  the  more  convenient 
relation  for  the  prediction  of  multi-phase,  multicomponent  equilibria: 


I  All  A  |  I  I 


A 

where  f^'s  are  the  fugacities  of  i  in  each  phase.  Restricting  our 
phases  to  liquid  and  vapor,  equation  (9)  would  result  to  the 
fundamental  criteria  for  vapor-liquid  equilibrium: 

qL  =  qv  (i=i,m)  do) 


The  general  vapor-liquid  equilibrium  problem  involves  a 
multicomponent  system  of  m  components  for  which  the  independent 
variables  are  T  and  P,  m-1  liquid  mole  fractions,  and  m-1  vapor 
mole  fractions.  Thus,  there  are  2m  independent  variables.  Application 


‘  .  * 
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of  the  phase  rule  shows  that  exactly  m  of  these  variables  must  be 
fixed;  while  the  remaining  m  variables  can  be  determined  by 
simultaneous  solution  of  m  equilibrium  relations  from  equation  (10). 
Specifying  either  the  liquid-phase  compositions  or  the  vapor-phase 
compositions  and  either  T  or  P  would  amount  to  fixing  m  independent 
variables.  Hence,  the  remaining  m  variables  are  subject  to 
calculation. 

Fugacity  is  expressed  for  each  phase  as: 

Liquid:  x^  P  (11) 

Vapor:  f.V  =  <J>.V  y.  P  (12) 

At  equilibrium: 

xi  P  =  yi  P  (13) 

This  introduces  the  composition  x.  and  y.  into  the  equilibrium  relation, 

A 

but  neither  are  explicit  because  the  f^'s  are  functions  of  composition. 

A 

In  addition,  <j>. 's  are  functions  of  T  and  P.  Thus,  equation  (13) 
represents  m  relationships  connecting  T,  P,  x.'s  and  y. 's.  For 

A 

computer  solution  of  this  equation  it  is  essential  that  the  <|>.  's  be 
expressed  analytically  as  functions  of  T,  P,  and  composition.  This 
requires  an  equation  of  state  which  accurately  represents  the 
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volumetric  properties  of  both  the  vapor  and  liquid  phases  throughout 
the  range  of  pressures,  temperatures ,  and  compositions  of  interest. 

4  23 

Most  of  the  computer  calculations  *  on  vapor-liquid 
equilibrium  use  the  two-parameter  Redl i ch-Kwong  equation  of  state 
in  determining  the  vapor  phase  fugacity  coefficient.  This  is 
probably  so  because  the  constants  of  the  equation  can  be  evaluated 
from  the  critical  temperature  and  pressure  for  the  pure  components 
constituting  the  system,  and  the  binary  interaction  coefficients 
k..  for  multicomponent  systems.  The  eight-parameter  Benedict-Webb- 

•  J 

,  .  24 

Rubin  (BWR)  equation  of  state  is  also  widely  employed,  but  a  large 
amount  of  volumetric  data  is  required  to  evaluate  the  constants. 

Such  constants  are  available  for  a  relatively  few  substances^  *  *  * 


Another  approach  in  solving  equation  (13)  is  by  using  the 
concept  of  activity  coefficients  in  the  liquid  phase.  Based  on  Raoult's 
law  and  the  concept  of  fugacity,  the  fugacity  of  a  component  in  a  liquid 
phase  is  given  by: 


A  I  0 

fi  =  Yi  xifi 


(14) 


where  y.  is  the  activity  coefficient  and  f^  the  standard  state  fugacity 

a  V 

coefficient  of  the  substance.  Retaining  the  same  relation  of  f.  , 
equation  (13)  will  become: 


Yi  Xi  fi°  -  *i  y,  p 


(15) 


. 

«  ■ 
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The  purpose  of  the  activity  coefficient,  y.  is  to  relate  the  unknown 
fugacity  to  its  standard  state  value,  f.®,  which  is  the  fugacity 
of  pure  i  at  the  temperature  and  pressure  of  the  system. 

The  activity  coefficient  is  related  to  the  Gibbs 
function  as  follows: 


GE 

rj*  “  I  (x-j  ^nYj ) 


(16) 


or, 


=  RT 


-  I 


X 


K 


3 (G  /RT) 


3X 


K 


T,  P,  x4 


(17) 


r  1  i ,  K 


Thus, for  a  given  T  and  P  and  composition,  y.  can  be  evaluated  if 

F 

G  is  known  as  a  function  of  composition. 

Activity  coefficients  can  also  be  calculated  from  experi¬ 
mental  vapor-liquid  equilibrium  data  by  the  use  of  the  following 
equation: 


Jtny.j  =  £n 


A 

+  £n  cj>.j 


£n  f . 


0 


(18) 


To  be  of  greater  use,  these  coefficients  are  often  correlated  by 
equations,  the  best  examples  being  the  van  Laar  equations  and  the 


Margules  equations. 


In  summing  up  the  analytical  solution  to  a  vapor-liquid 
equilibrium,  it  is  sufficient  to  say  that  given  a  suitable  equation 
of  state  with  known  constants  for  the  vapor  phase,  data  for  the 
standard  states,  and  an  appropriate  relation  of  activity  coefficients 
in  the  liquid  phase,  the  unknown  variables  can  be  determined. 


% 


13 


III  -  PREVIOUS  WORK 
A)  Carbonyl  Sulfide  (COS) 

The  existence  of  carbonyl  sulfide  was  not  confirmed  until 
29 

1867  when  Than  reported  it  by  the  reaction  of  carbon  monoxide 
and  sulfur  vapors.  Since  then,  numerous  investigators  have 
prepared  it  in  a  variety  of  ways.  A  comprehensive  review  of  these 

3 

methods  and  the  chemistry  involved  has  been  presented  by  Ferm  . 

(i)  Vapor  Pressure 

Earlier  studies  on  the  vapor  pressure  of  carbonyl 

30  31  32 

sulfide  were  made  by  Ilosvay  ,  Hempel  ,  and  Stock  and  Kuss  . 

These  investigations  however,  were  conducted  only  for  a  limited 

33 

number  of  temperatures.  The  International  Critical  Tables  report 

the  work  of  Ilosvay  but  some  contradictions  in  the  reported  data 

34 

were  observed.  In  1936,  Kemp  and  Giauque  measured  COS  vapor 
pressures  from  161.8  K  to  223.8  K.  The  data  were  fitted  to  the 
fol lowing  equation : 

log1Q(P)  =  -(1318.260/T)  +  10.15309  -  0.0147784T 

+  0.00001 8838T2  (19) 

where  P  is  the  absolute  pressure  in  cm.  Hg.,  and  T  is  the  absolute 

35 

temperature  in  K.  Later,  Honig  and  Hook  developed  another  equation 


' 


. 
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to  fit  these  vapor  pressure  observations.  The  Matheson  Gas  Products 

36 

Company  in  their  latest  publication  on  COS  properties  presented  a 
vapor  pressure  diagram  based  on  the  equation  of  Kemp  and  Giauque,  but 
extended  its  use  to  temperatures  well  beyond  the  range  for  which  it 
was  developed. 

(ii)  Density 

The  density  of  saturated  liquid  COS  at  -87°C  was  measured 

32 

to  be  1.24  g/cc  by  Stock  and  Kuss  as  compared  to  a  value  of  1.30  g/cc 

37 

suggested  by  Pearson,  Robinson,  and  Trotter  at  the  same  conditions. 
Orthobar ic  densities  for  liquid  COS  were  obtained  by  Partington  and 

oo 

Neville  from  -99.5°C  to  100. 9°C,  however,  their  data  showed  a 

°2 

disagreement  of  up  to  4%  when  compared  to  the  work  of  Stock  and  Kuss° 

37 

and  Pearson  et  al  .  The  only  vapor  density  measurement  for  COS 

34 

was  reported  by  Kemp  and  Giauque  .  At  25°C  and  atmospheric  pressure, 
the  value  was  reported  to  be  2.4849  g/cc. 

(iii)  Critical  Properties 

The  critical  temperature  of  COS  was  reported  to  be  105°C 
30  31  • 

by  Ilosvay  and  Hempel  ,  but  it  was  measured  to  be  102. 2°C  by 

38 

Partington  and  Neville  .  From  the  vapor  pressure  data  of  Stock  and 

39 

Kuss,  Kobe  and  Lynn  estimated  the  critical  pressure  to  be  61  atm. 

A  critical  pressure  estimation  made  by  Lydersen^,  however,  indicated 


a  value  of  65  atm. 


-  • 
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(iv)  Refractive  Index 

The  refractive  index  of  liquid  COS  was  measured  by  Francis^ 
42 

and  Wong  and  Anderson  ;  the  latter  reported  the  molar  refractivity 
to  be  11.29  cc/g. -mole. 

(v)  Other  Properties 

43 

The  phase  diagram  for  COS  was  presented  by  Babb  ,  and 

the  solubility  of  COS  with  some  organic  solvents  were  measured  by 

.  44 

Kir  yanova  and  Pinsker  . 

In  addition  to  vapor  pressure  and  vapor  density  measure- 
34 

ments,  Kemp  and  Giauque  also  determined  its  melting  point,  boiling 
point,  heat  capacities,  and  heats  of  fusion  and  vaporization.  Other 

thermodynamic  properties  for  COS  have  been  calculated  by  Kobe  and 

39  45  46 

Lynn  ,  Kobe  and  Long  ,  McBride  and  Gordon  ,  Zandler,  Watson,  and 

47  48 

Eyring  ,  and  Cross  . 

Because  COS  occurs  as  a  contaminant  in  many  refinery 

gases,  its  corrosive  properties  are  of  interest.  Nelson,  Veal,  and 
49 

Heinrich  verified  that  COS  does  not  produce  a  tarnished  strip 

50 

in  their  copper  corrosion  test.  However,  Freise  discovered  that 
COS  is  corrosive  toward  concrete.  Except  for  these  findings,  in¬ 
formation  is  not  available  on  the  action  of  COS  towards  other  metals. 

B)  Propane 


Information  on  the  critical  properties,  vapor  pressure 


•  '  •  •  •  •  i  .  _ 
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characteristics,  density  and  other  thermodynamic  properties  of 

51  -56 

propane  have  been  reported  by  many  authors  .  Its  behaviour  with 

other  hydrocarbons  and  nonhydrocarbons  is  extensively  reported  in 
57-75 

the  literature  .  A  detailed  review  on  this  work  will  not  be 
included  here  for  it  is  beyond  the  scope  of  this  study. 

C)  Carbonyl  Sulfide-Propane  Binary 

No  information  is  available  on  this  system.  And  as 
mentioned  earlier,  no  studies  have  been  conducted  on  the  vapor- 
liquid  equilibrium  of  carbonyl  sulfide  with  other  substances. 


. 
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IV  -  THE  EXPERIMENTAL  STUDY 

The  objective  of  the  experimental  program  was  to  obtain 
pressure-composition  diagrams  for  the  carbonyl  sulfide-propane 
system  simultaneously  with  the  co-existing  phase  densities. 

The  densities  of  the  co-existing  phases  were  calculated 
from  refractive  index  measurements  which  were  taken  twice  for  each 
phase.  Vapor  and  liquid  phase  compositions  were  analyzed  using  a 
gas  chromatograph.  Two  samples  of  each  phase  were  analyzed  three 
times  and  the  average  of  the  six  analyses  was  accepted. 

A)  Experimental  Apparatus 

A  detailed  description  of  the  equipment  design  has  been 

7fi  77 

presented  in  earlier  publications  *  .  Although  no  significant 

modification  was  made  to  the  original  equipment,  it  is  considered 
desirable  to  mention  some  of  the  more  pertinent  features  here. 

A  schematic  diagram  of  the  apparatus  and  its  associated 
equipment  is  shown  in  Figure  2  and  the  essential  features  of  the 
cell  design  are  shown  in  Figure  3.  The  cell,  machined  from  4  inch 
diameter  type  316  stainless  steel,  consists  of  three  parts,  two 
cylinder  piston  end  sections  and  a  central  windowed  section.  The 
three  sections  are  bolted  together  with  high  pressure  seals  between 
them  made  by  flattened  Teflon  0-rings.  On  assembly,  it  has  an  over- 
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FIG.  2  SCHEMATIC  DIAGRAM  OF  THE  CELL  AND  ITS  ASSOCIATED  EQUIPMENT 


19 


FIG.  3  ESSENTIAL  FEATURES  OF  THE  EQUILIBRIUM  CELL  DESIGN 


20 


all  length  of  19  inches.  Each  piston  has  a  4  inch  travel  and  is 
confined  to  its  respective  cylinder.  A  low  viscosity  silicone  oil 
is  employed  as  the  hydraulic  fluid.  This  provides  a  means  for 
varying  the  cell  volume  and  pressure.  The  piston  seal  between  the 
oil  chamber  and  the  central  windowed  section  is  effected  by  four 
neoprene-O-rings.  The  original  design  made  use  of  a  316  type 
stainless  steel  piston.  However,  later  indications  showed  that 
wearing  out  of  the  cell's  inner  surface  resulted  from  metal  to  metal 
contact  during  the  travelling  of  the  pistons,  especially  when 
working  at  high  pressure.  To  resolve  the  problem,  the  pistons  were 
re-machined  from  Teflon. 

The  horizontal  section  A-A  through  the  center  of  the  cell 
shows  how  the  window,  a  pyrex  disk  1-1/8  inches  in  diameter  and  1  inch 

_5 

thick  with  its  faces  flat  and  parallel  to  ±4  x  10  inches,  and  a 
type  316  stainless  steel  mirror  form  the  boundaries  of  a  30°  prism. 

The  window  and  mirror  are  held  in  place  by  bolted  cover  plates  and 
sealed  against  lapped  surface  with  1/16  inch  thick  glass-filled 
Teflon  gaskets. 

The  horizontal  section  B-B  shows  one  of  the  sampling 
valves,  the  pressure  transducer  port,  and  the  circulation  line 
outlet.  The  sampling  valves  are  internally  mounted  Circle  Seal  MV-92 
micro-metering  valves  with  provision  made  for  flushing  out  and 
evacuation  of  the  low  pressure  sampling  line.  The  circulation 
system  used  to  attain  equilibrium  consists  of  an  internal  check 
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valve  in  the  central  section  and  an  external  1/8  diameter  stainless 
steel  line  with  a  shut-off  valve. 

Cell  Temperature  Control 

The  temperature  of  the  cell  is  maintained  by  two  7  inch 
diameter  aluminum  environmental  control  shrouds  which  are  placed 
over  the  end  of  the  cell  as  shown  in  Figure  2.  Each  shroud  includes 
four  150  watt  strip  heaters  and  a  20  foot  cooling  coil.  A  Thermae 
proportional  plus  integral  controller  operates  the  heaters  in  both 
shrouds  and  controls  the  cell  temperature  to  within  ±  0.1°F  of  the 
setpoint. 

When  working  at  isotherms  below  ambient  temperature,  the 
cell  is  cooled  by  circulating  ethylene  glycol-water  mixture 
refrigerated  by  freon-12. 

The  cell  is  capable  of  operating  from  2°F  to  250°F  in 
normal  ambient  room  conditions. 

Optical  System 

A  telescope  (Gaertner  M523  aperture  28  mm.,  f.l.  250  mm.) 
with  an  Abbe-Lamont  autocol 1 imating  eyepiece  (Gaertner  1372)  is 
mounted  on  a  10  inch  diameter  precision  rotary  table  (Karl  Kneise 
RTHP-10)  with  a  two  second  vernier  handwheel  and  a  maximum  table 
error  of  10  seconds  in  360  degree  of  rotation.  The  table  is  centered 
under  a  point  on  the  front  surface  of  the  pyrex  cell  window.  Mono- 

O 

chromatic  light  at  a  wavelength  of  6328  A  is  provided  by  a  Helium- 
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Neon  gas  laser  and  is  transmitted  to  the  eyepiece  by  a  glass  fiber 
light  guide. 

The  refraction  angle  y  shown  in  Figure  1  is  measured  by 
autocol 1 imation .  The  measured  angles  have  a  repeatability  of 
±4  seconds  and  a  table  error  of  ±2  seconds  over  the  measured  range. 

_5 

The  refractometer  has  an  over-all  possible  error  of  ±6  x  10  in 
the  refractive  index  and  a  possible  range  of  1.0  to  1.7. 

Operational  Procedure 

The  conventional  way  of  working  with  binary  systems 
involves  charging  the  heavier  component  first  following  by  adding  the 
lighter  component.  In  this  study,  the  carbonyl  sulfide  (the  lighter 
component)used  has  a  maximum  purity  of  97.6  mole  percent  in  the  liquid 
phase  and  as  such  had  to  be  purified  before  any  binary  data  could  be 
taken.  The  bulk  of  the  impurities  consisted  of  lighter  substances 
which  could  be  effectively  removed  by  fractional  distillation.  For 
convenience,  the  carbonyl  sulfide  was  charged  first  and  was  fractionally 
distilled  inside  the  cell.  Propane  was  subsequently  added  to  make 
up  the  pressure  to  the  desired  level.  This  was  achieved  by  boosting  the 
feeding  pressure  of  propane  in  the  cylinder  with  a  hot  water  bath.  A 
propane  cylinder  pressure  of  at  least  65  psia  greater  than  the  cell 
pressure  was  considered  adequate  to  prevent  contamination  of  the  cylinder 
with  the  contents  of  the  cell. 

The  two  pistons  in  the  end  sections  were  hydraulically 
driven  by  a  motorized  dual  cylinder  high  pressure  pump  with  an 
opposed  transmission  (Ruska  Model  2248  WII).  This  allowed  both 
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pistons  to  be  moved  simultaneously  up  or  down  maintaining  the 
volume  between  them  constant.  The  possible  working  volume  ranges 
from  10  cc  to  174  cc. 

Mixing  and  equilibration  were  attained  by  running  the 
pistons  up  and  down  simultaneously.  On  the  upstroke,  the  internal 
check  valve  closes  and  the  fluid  is  forced  through  the  external 
circulation  line  and  sprayed  into  the  upper  cylinder  through  a 
nozzle.  On  the  downstroke,  the  circulation  line  valve  is  closed 
and  the  fluid  passes  through  the  internal  check  valve.  Five  to 
ten  cycles  were  usually  sufficient  to  attain  equilibrium. 

The  point  at  which  equilibrium  is  achieved  can  be  verified 
by  observing  the  line  images  in  both  phases  through  the  auto- 
collimating  telescope.  These  images  are  very  sensitive  to  composition 
and  temperature  fluctuations  and  become  stable  when  the  equilibrium 
state  is  attained.  The  angle  y  was  measured  by  observing  the 
location  of  the  line  images  obtained  through  normal  light  incidence 
on  the  glass  window  and  the  mirror. 

The  temperature  was  measured  with  an  i ron-constantan 
thermocouple  sheathed  in  316  stainless  steel  with  its  tip  exposed  to 
the  fluid  in  the  prism.  The  pressure  was  measured  with  a  3000  psi 
316  stainless  steel  bourdon  tube  Heise  Gauge  or  a  1500  psi  strain 
gauge  transducer  (Consolidated  Electronics). 

Samples  of  liquid  or  vapor  were  expanded  to  a  pressure  of 
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0.2  atmospheres  through  a  micro-metering  valve  into  the  evacuated 
line  connected  to  the  gas  chromatograph  gas  sampling  valve.  The 
pressure  in  the  sampling  line  was  measured  with  a  differential 
pressure  transducer.  Two  samples  of  each  phase  were  taken  and 
triplicate  chromatographs  were  run  on  each  sample.  The  average  size 
of  each  sample  corresponds  to  a  0.2  percent  depletion  of  an  average 
load.  The  pressure  drop  associated  with  the  taking  of  a  sample  was 
always  less  than  3  psia.  For  each  point  the  six  analyses,  three 
from  each  of  two  samples,  were  generally  repeatable  to  within  0.2 
mole  percent  and  their  accuracy  is  probably  0.3  percent.  The 
pressure  is  believed  known  to  ±3  psia  and  the  temperature  to  ±0.1 °F. 


The  gas  chromatograph  used  was  a  Hewlett-Packard  Model 


700  machine  equipped  with  a  thermal  conductivity  cell  detector. 

The  signal  from  the  gas  chromatograph  was  detected  by  an  on-line 
IBM-1800  Computer  and  peak  integrations  were  obtained  by  using 
available  software.  The  results  of  the  analysis  were  printed  out  on 
a  teletype. 

Operating  conditions  for  the  gas  chromatograph  were  as 

follows: 


Column 


Porapak  QS 


Column  Size 


1/8"  dia.  x  6 


Carrier  Gas 


Helium  at  20  c.c./min. 
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Column  Temperature 
Filament  Current 
Detector  Temperature 
Area  Integration 


75°C 
180  ma. 
200°C 

Electroni c 


The  conditions  given  above  were  experimental ly  determined 
to  be  optimum  considering  separation  and  time  of  analysis.  The  gas 
chromatograph  was  calibrated  using  the  pure  components  only  and  the 
response  was  found  to  be  linear  for  the  range  of  the  sample  loop 
pressures  encountered.  The  resulting  calibration  curve  is  shown  in 
Figure  20.  The  response  factor  of  propane  relative  to  carbonyl 
sulfide  was  calculated  to  be  0.995. 


B)  Materials 

Matheson  instrument  grade  carbonyl  sulfide  was  used. 
Chromatographic  analysis  indicated  it  to  contain  97.6  mole  percent 
COS  in  the  liquid  phase  and  84.1  mole  percent  COS  in  the  gas  phase, 
the  major  impurity  being  carbon  dioxide.  Because  of  the  impurity, 
further  purification  was  necessary,  and  it  was  carried  out  in  two 
stages.  Liquid  COS  was  drawn  into  a  500  ml.  stainless  steel  sample 
bomb  (Hoke  DS  5-52),  connected  to  a  3000  psi  stainless  steel  Heise 
gauge.  After  isolation  from  the  main  cylinder  the  contents  of  the 
sample  bomb  were  cooled  by  liquid  nitrogen.  When  the  gauge 
registered  a  pressure  of  about  10  psi,  the  supply  of  nitrogen  was  shut 
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off,  and  the  sample  was  fractionally  distilled  by  venting  the  vapor 
phase.  To  minimize  the  loss  of  COS,  the  fractionation  was  stopped 
when  the  pressure  reached  about  120  psig.  This  step  was  carried 
out  four  to  five  times.  Analysis  of  the  vapor  phase  after  the  final 
fractionation  indicated  it  to  contain  approximately  98.5  mole  percent 
COS.  Further  fractionation  made  very  small  improvement  in  the 
purity. 

The  second  stage  in  the  purification  was  done  inside  the 
equilibrium  cell.  A  liquid  sample  of  about  120  cc,  purified  from 
the  first  stage,  was  charged  to  the  cell.  Since  the  cell  was  equipped 
with  a  micro-metering  valve,  the  vapor  phase  could  be  bled  off  slowly. 

The  pistons  in  the  end  sections  were  adjusted  occasionally  in  such  a 
manner  that  the  interface  of  the  liquid  and  vapor  was  always  visible 
in  the  window,  and  vaporization  of  the  liquid  could  be  kept  to  a 
minimum.  When  about  40  cc  of  liquid  were  left,  the  vapor  phase  was 
analyzed.  It  was  found  to  contain  about  99.5  mole  percent  COS. 

Subsequent  bleeding  off  had  a  negligible  effect  on  the  purity.  The 
COS  used  in  the  study  was  thus  of  a  99.5  percent  minimum  purity.  The 
second  stage  purification  generally  took  about  4  to  5  hours. 

The  propane  used  was  Matheson  Instrument  grade.-  It  was 
analyzed  and  found  to  contain  better  than  99.9  mole  percent  in  the 
liquid  phase  and  98.0  mole  percent  in  the  vapor  phase,  the  major  impurities 
being  methane  and  nitrogen.  Due  to  this  purity  the  liquid  phase  was 
used  in  the  study. 


. 
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V  -  EXPERIMENTAL  RESULTS 

The  phase  behaviour  of  the  carbonyl  sulfide-propane 
binary  system  was  studied  at  20°,  60°,  100°,  140°,  and  180°F, 
with  pressures  ranging  from  the  vapor  pressure  of  propane  to  the 
vapor  pressure  of  carbonyl  sulfide.  An  average  of  10  points, 
generally  equally  spaced  in  pressure  between  pure  propane  and  pure 
COS  vapor  pressures  were  taken  for  each  isotherm. 

The  equilibrium  phase  compositions  for  each  of  the  five 
isotherms  were  analyzed.  The  refractive  indices  of  the  co-existing 
phases  were  measured  and  the  molar  volumes  calculated  using  a  molar 
averaged  Lorentz-Lorenz  refractivity.  The  refractive  indices  of 
pure  propane  and  pure  COS  were  also  measured  for  each  of  the  isotherms 
from  saturation  pressure  up  to  about  1000  psia.  For  COS,  the 
refractive  index  measurement  was  extended  to  230°F.  This  provided 
sufficient  data  to  plot  saturation  pressure  versus  refractive  index, 
from  which  an  estimate  of  the  critical  pressure  could  be  obtained. 

The  vapor  pressure  curve  for  COS  was  plotted  from  20°F  to  the  estimated 
critical  point. 

A)  Refractive  Indices  of  Carbonyl  Sulfide  and  Propane 

Experimentally  determined  refractive  indices  of  the 
compressed  liquid  for  pure  COS  and  propane  are  plotted  with  pressure. 
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and  are  presented  in  Figures  4  and  5,  respectively.  These  data 
are  also  tabulated  in  Tables  1  and  2  in  the  Appendix. 

B)  Refractive  Indices  of  Carbonyl  Sulfide  at  Saturated  Conditions 

Table  3  in  the  Appendix  shows  the  experimental ly  measured 
refractive  indices  of  saturated  carbonyl  sulfide  from  20°F  to  230°F. 
Figure  6  gives  a  plot  of  the  saturation  pressure  versus  the  re¬ 
fractive  index  for  the  vapor  and  the  liquid. 

The  critical  point  estimated  from  Figure  6  is  plotted  along 
with  the  vapor  pressure  data  obtained  for  carbonyl  sulfide  in 
Figure  7. 

C)  Co-existing  Phase  Properties  of  the  Carbonyl  Sulfide-Propane  System 

Figures  8,  9,  10,  11,  and  12  show  the  plots  of  the 
experimental ly  determined  phase  compositions  with  pressure.  The 
data  indicate  the  distinct  S-shaped  configuration  of  the  P-x  diagrams. 

Plots  of  the  measured  refractive  indices  are  shown  in 

Figures  13  and  14.  The  molar  -refractivi ties  of  propane  were  calculated 

using  the  volumetric  data  predicted  by  the  BWR  equation  with  improved 
28 

coefficients  ,  and  the  measured  refractive  indices.  For  COS,  the 
necessary  volumetric  data  were  interpolated  from  those  values  reported 
by  Partington  and  Neville  .  The  molar  volumes  of  the  co-existing 
phases  were  calculated  using  a  molar  average  Lorentz-Lorenz  re- 
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FIG.  4  EFFECT  OF  PRESSURE  AND  TEMPERATURE  ON  THE 
REFRACTIVE  INDEX  OF  CARBONYL  SULFIDE 
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FIG.  5  EFFECT  OF  PRESSURE  AND  TEMPERATURE  ON  THE 
REFRACTIVE  INDEX  OF  PROPANE 
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FIG.  6  PRESSURE  VS.  REFRACTIVE  INDEX  PLOT  OF  CARBONYL  SULFIDE 
AT  SATURATED  CONDITIONS  FOR  CRITICAL  POINT  DETERMINATION 
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FIG.  7  VAPOR  PRESSURE  CURVE  FOR  CARBONYL  SULFIDE 
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.  3  PRESSURE-EQUILIBRIUM  PHASE  COMPOSITION  DIAGRAM  FOR  THE 
CARBONYL  SULFIDE-PROPANE  SYSTEM  AT  20°F 
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9  PRESSURE-EQUILIBRIUM  PHASE  COMPOSITION  DIAGRAM  FOR  THE 
CARBONYL  SULFIDE-PROPANE  SYSTEM  AT  60°F 
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FIG.  10  PRESSURE-EQUILIBRIUM  PHASE  COMPOSITION  DIAGRAM  FOR  THE 
CARBONYL  SULFIDE-PROPANE  SYSTEM  AT  100°F 
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11  PRESSURE-EQUILIBRIUM  PHASE  COMPOSITION  DIAGRAM  FOR  THE 
CARBONYL  SULFIDE-PROPANE  SYSTEM  AT  140°F 
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FIG.  12  PRESSURE-EQUILIBRIUM  PHASE  COMPOSITION  DIAGRAM  FOR  THE 
CARBONYL  SULFIDE-PROPANE  SYSTEM  AT  180°F 
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FIG.  13  PRESSURE-EQUILIBRIUM  PHASE  REFRACTIVE  INDEX  DIAGRAM 

FOR  THE  CARBONYL  SULFIDE-PROPANE  SYSTEM  AT  140°F  AND  180°F 
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FIG.  14  PRESSURE-EQUILIBRIUM  PHASE  REFRACTIVE  INDEX  DIAGRAM  FOR 

THE  CARBONYL  SULFIDE-PROPANE  SYSTEM  AT  20°F,  60°F  AND  100°F 


f r act i v i t y : 
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These  are  shown  in  Figures  15  and  16. 

K-factor  versus  pressure  curves  are  presented  in 
Figures  17  and  18. 


The  above  data  are  also  tabulated  in  the  Appendix. 
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FIG.  15  PRESSURE-EQUILIBRIUM  PHASE  MOLAR  VOLUME  DIAGRAM  FOR  THE 
CARBONYL  SULFIDE-PROPANE  SYSTEM  AT  140°F  AND  180°F 
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FIG.  16  PRESSURE-EQUILIBRIUM  PHASE  MOLAR  VOLUME  DIAGRAM  FOR  THE 
CARBONYL  SULFIDE-PROPANE  SYSTEM  AT  20°F,  60°F  AND  100°F 
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FIG.  17  EQUILIBRIUM  RATIOS  FOR  CARBONYL  SULFIDE-PROPANE  SYSTEM  AT 
140°F  AND  1 80°F 
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FIG.  18  EQUILIBRIUM  RATIOS  FOR  CARBONYL  SULFIDE-PROPANE  SYSTEM  AT 
20°F»  60°F,  AND  100°F 


- 


45 


VI  -  DISCUSSION  OF  RESULTS 

A)  The  System 

In  view  of  the  fact  that  very  little  is  known  about  the 
nature  of  carbonyl  sulfide  some  comments  about  its  physical  and 
chemical  behaviour  with  the  experimental  equipment  seems  relevant. 

Carbonyl  sulfide  has  been  described  as  a  colorless  gas 
with  no  information  available  on  the  color  of  the  liquid.  However, 
during  this  investigation  the  liquid  was  observed  to  be  bright  red 
in  color.  One  might  speculate  about  the  possible  cause  of  this 
color  wherein  at  least  two  plausible  postulates  could  be  examined. 

As  a  first  possibility  liquid  COS  may  in  fact  be  red  or  alternati vely 
the  color  maybe  due  to  some  impurity  generated  by  a  chemical  reaction 
of  COS  with  the  materials  of  construction  of  the  equilibrium  cell. 

Two  possible  reactions  of  COS  with  iron  (Fe)  may  be  written  down  as: 

COS  +  Fe  -*  FeS  +  COt  (21) 

COS  +  Fe  +  FeCO  +  S*  (22) 

Equation  (21)  suggests  evolution  of  carbon  monoxide  but  no  evidence 
of  this  was  observed  during  vapor  phase  analysis  on  the  gas  chromatograph. 
Equation  (22)  appears  to  be  the  more  plausible  explanation  since  FeCO 
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is  known  to  be  red  and  a  black  solid  deposit  (possibly  impure  sulfur) 
was  noticed  on  the  inner  surface  of  the  equilibrium  cell.  More 
analytical  work  is  however  needed  to  confirm  this  hypothesis. 

The  red  color  mentioned  above  posed  considerable 
difficulties  in  measuring  the  refractive  index  of  the  liquid  phase 

O 

mainly  due  to  the  fact  that  the  laser  light  source  used  (6328  A)  was 
red  in  color  and  the  ray  of  light  transmitted  through  the  liquid  is 
to  a  great  extent  absorbed. 

B)  Determination  of  the  Critical  Point  for  Carbonyl  Sulfide 

The  critical  point  was  determined  by  plotting  vapor 
pressures  versus  refractive  indices  of  the  saturated  liquid  and  vapor 
phases.  The  critical  temperature  of  COS  was  estimated  by  reading  off 
the  temperature  corresponding  to  the  critical  pressure  from  a  smoothed 
plot  of  the  experimentally  measured  vapor  pressure  data  versus  temp¬ 
erature  (Figure  7).  The  critical  pressure  and  temperature  obtained 
correspond  to  960  psia  and  235°F  respectively. 

Figure  19  shows  a  plot  of  the  reduced  pressure  Pr  versus 

the  reduced  temperature  T^  of  saturated  COS.  Experimental  data  are 

compared  to  a  similar  plot  obtained  from  generalized  tables  of 

saturated  fugacity  coefficients  as  a  function  of  Pr  and  Tr  at  different 
79 

values  of  Zc  .  Excellent  agreement  is  obtained  between  the 

generalized  plot  of  Z  =0.29  and  the  experimental  values  as  indicated 
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REDUCED  TEMPERATURE 


FIG.  19  REDUCED  PRESSURE  VS.  REDUCED  TEMPERATURE  PLOT  AT 
SATURATED  CONDITIONS  FOR  CARBONYL  SULFIDE 
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in  Figure  19.  The  critical  compresibi 1 ity  factor  Z  was  calculated 

c 

on 

to  be  0.287  from  the  Curl  and  Pitzer's  correlation: 

Zc  =  1.28o)  +  3.41 

The  critical  properties  estimated  from  this  work  indicate 

TO  T1  T8  T9 

some  discrepancies  when  compared  to  the  work  of  other  investigators 

The  difference  may  be  attributed  to  the  purity  of  the  COS  used  in 
each  case  and  the  techniques  employed  to  obtain  the  critical  properties. 

C)  Equilibrium  Phase  Properties  of  the  COS-Cq  System 

The  equilibrium  phase  composition  data  obtained,  as 
shown  on  pressure-composition  diagrams,  suggest  that  no  azeotropes 
are  formed  for  this  system  in  the  temperature  range  investigated. 

Figures  8,  9,  10,  11,  and  12  show  an  elongated  S-shape  phase  diagram. 

An  examination  of  these  figures  indicates  possibility  of  azeotropic 
behavior  at  the  carbonyl  sulfide  rich  end,  the  trend  becoming  more 
pronounced  at  higher  temperatures. 

The  scattered  points  on  the  plots  of  pressure  versus 
refractive  index  and  molar  volume  shown  in  Figures  13,14  and  15,16 
are  indicative  of  the  difficulties  encountered  in  obtaining  the 
refractive  index  readings. 

Calculated  equilibrium  ratios  based  on  the  experimental 
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vapor-liquid  equilibrium  compositions  are  of  a  conventional  shape 
when  plotted  with  pressure.  All  of  the  isotherms  exhibit  a  minimum 
value  of  K  for  the  heavier  component.  These  are  shown  in  Figures 


17  and  18. 
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VII  -  PREDICTION  OF  THE  DATA  AND  THERMODYNAMIC 

CONSISTENCY  TEST 


General 

A  number  of  predictive  methods  are  available  for 
calculating  multicomponent  vapor-liquid  equilibrium.  Three  of  the 
more  widely  used  are  those  of  a)  Chueh  and  Prausnitz^  b)  Benedict, 
Webb,  and  Rubin  and  c)  Soave-Redl ich-Kwong  .  Each  of  these  is 
available  in  the  form  of  a  computer  program.  Basically  each  method 
utilizes  the  concept  that  f^  =  f.^  at  equilibrium.  At  a  fixed 
temperature  and  liquid  composition,  the  pressure  and  the  vapor 
composition  that  satisfies  the  above  condition  are  calculated.  The 
Chueh-Prausni tz  method  was  used  exclusively  for  this  investigation. 

The  modified  version  of  the  Chueh-Prausni tz  phase 
behavior  program  called  "CHUBIN"  and  the  program  called  "FITTING" 
which  estimates  the  value  of  the  interaction  constant  a-^  and  the 
dilation  constant  n  were  used. 

Chueh-Prausnitz  Parameters 

Owing  to  lack  of  binary  data  involving  carbonyl  sulfide  with 
hydrocarbons,  no  information  is  available  for  a  number  of  parameters 
needed  in  the  "CHUBIN"  and  "FITTING"  programs  to  predict  the  vapor- 
liquid  equilibrium.  Some  of  these  parameters  were  estimated  using  the 
critical  properties  of  the  pure  components  and  some  were  determined 
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by  using  the  data  obtained  at  100°F.  The  constants  were  estimated 
so  that  the  100°F  isotherm  could  be  predicted.  These  constants 
were  subsequently  used  to  predict  the  phase  behavior  at  the  other 
temperatures.  The  estimated  parameters  are  listed  in  Table  5. 

Prediction  of  the  Critical  Properties  of  the  Binary 

The  critical  temperature,  pressure,  and  molar  volume  of 
the  binary,  with  mole  fraction  compositions  incremented  from  0.0  by 
0.05  to  1.00  were  predicted  by  the  Chueh-Prausni tz  computer  program. 

The  predictions  are  tabulated  in  Table  7.  The  critical 
locus  exhibited  a  minimum  temperature  of  194. 1°F  with  a  corresponding 
pressure  of  559  psia. 

Thermodynamic  Consistency  Test 

The  experimental  data  were  tested  for  thermodynamic 

consistency  by  using  a  method  presented  by  Chueh,  Muirbrook,  and 
78 

Prausnitz  .  The  test  was  derived  from  the  Gibbs-Duhem  equation 
which  for  a  binary  system  at  constant  temperature  is: 

x-|d£nf-j  +  =  v^dP/RT 

By  appropriate  thermodynamic  substitution  and  mathematical  man¬ 
ipulation,  the  authors  show  that  a  point  by  point  consistency  test 
can  be  arrived  at  as  follows: 
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Area  I  +  Area  II  +  Area  III  = 


(f)-.  P  <f>o  Kp 

[£n  K-j  +  £n  — g — ^  +  x2(£n  +  Jin  )  ]  at  Xp 
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The  three  areas  were  numerically  integrated  by  using  a 

pi 

method  of  assigned  sample  points  at  unequal  intervals  .  The 
fugacity  coefficients,  <j>.  's  were  calculated  from  the  Chueh- 
Prausnitz  program,  and  the  equilibrium  ratios,  K. 's  were  estimated 
from  the  experimental  phase  compositions. 


Area 


111  =w 


Due  to  possible  uncertainties  in  the  evaluation  of  the 

fugacity  coefficients  through  the  use  of  an  equation  of  state, 

78 

the  authors  judged  that  five  per  cent  inconsistency  constitutes 
a  very  good  thermodynamic  consistency  test  of  the  data. 

Table  8  shows  a  comparison  of  the  three  areas,  the  left 
hand  side  and  right  hand  side  of  the  equation,  and  the  percentage 
inconsistency  for  the  five  isotherms  all  calculated  at  Xp=1.0. 

The  test  indicates  an  increasing  inconsistency  at  elevated 
temperatures.  This  might  be  due  to  three  possible  reasons.  Firstly, 
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the  technique  used  in  evaluating  fugacity  coefficients  at  high 
temperatures,  especially  near  the  critical  region  might  contribute 
a  significant  error.  Secondly,  the  possibility  of  large  round  off 
errors  introduced  by  the  nature  of  fugacity  coefficients  and  the 
equilibrium  ratios  for  each  component  are  very  close  to  each  other 
which  suggests  that  the  ratio  would  be  always  close  to  unity.  The 
natural  logarithm  of  these  ratios  would  be  close  to  zero  where  small 
errors  in  the  calculation  would  get  magnified.  Finally,  the  data 
itself  could  have  some  inconsistency. 


• 
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VIII  -  SUMMARY  AND  CONCLUSIONS 

The  results  and  conclusions  of  the  work  undertaken  and 
reported  herein  may  be  summarized  as  follows: 

1.  The  equilibrium  phase  compositions  and  refractive 
indices  have  been  measured  for  the  carbonyl  sulfide-propane  binary 
system  at  20°,  60°,  100° ,  140°,  and  180°F  and  at  pressures  between 
the  vapor  pressure  of  propane  and  carbonyl  sulfide. 

2.  The  refractive  indices  of  pure  carbonyl  sulfide 

and  pure  propane  were  measured  at  each  of  the  above  temperatures  to 
pressures  of  about  1000  psia. 

3.  The  refractive  indices  of  saturated  liquid  and  saturated 
vapor  carbonyl  sulfide  were  measured  at  11  temperatures  between  20° 
and  230°F. 

4.  The  vapor  pressure  of  carbonyl  sulfide  has  been  measured 
from  20°F  to  230°F. 

5.  The  refractive  index  measurements  and  the  vapor  pressure 
and  temperature  measurements  of  the  coexisting  vapor  and  liquid 
phases  for  pure  carbonyl  sulfide  have  been  used  to  estimate  the 
following  properties  for  carbonyl  sulfide: 

Critical  Pressure  960  ±  5  psia 

Critical  Temperature  235  ±  1°F 


Acentric  Factor 


0.056 


. 
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The  critical  compressibility  factor  was  calculated 

3 

to  be  0.287  and  the  critical  volume  was  calculated  to  be  2.227  ft  /lb-mole. 

5.  Numerical  values  have  been  obtained  for  the  binary 
parameters  required  in  using  the  Chueh-Prausnitz  correlation  for 
carbonyl  sulfide-propane  binaries,  including  k. the  deviation 

'  sJ 

from  the  geometric  mean  for  T  ;  xi •  the  parameter  characteristic 

c  i  j  i  j 

of  the  i-j  interaction  for  T  ;  v..,  the  parameter  characteristic 

c  1  j 

of  the  i-j  interaction  for  V  ;  a..,  the  interaction  parameter 

C  1 J 

characteristic  of  the  i-j  interaction  in  the  liquid  solution 
at  each  of  the  experimental  temperatures;  and  and  ^  for  both 
liquid  and  vapor  required  in  the  Redl ich-Kwong  equation  of  state. 
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NOMENCLATURE 

Angstrom 

Centigrade 

Fahrenheit 

Gibbs  free  energy 

Vapor-liquid  equilibrium  constant 

Molecular  weight 

Pressure 

Universal  gas  constant 

Lorentz-Lorenz  refractivity 

Entropy 

Temperature 

Internal  energy 

Molar  volume 

Compressibility  factor 

Fugacity 

Refractive  index 

Mole  fraction  in  the  liquid 

Mole  fraction  in  the  vapor 


Greek  letters 


Redl ich-Kwong  constants 
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a 

3 

B 

Y 

Y 
6 

P 

n 

T 

Y 

v 

03 


van  Laar  interaction  parameter 
Angle  refraction  of  the  window 

Correction  for  critical  properties  in  the  Chueh-Prausnitz 
computer  program 
Activity  coefficient 

Minimum  deviation  angle  measured  by  autocol 1 imation 

Prism  angle 

Density 

Fugacity  coefficient 

Dilation  constant  in  the  Chueh-Prausnitz  program 
Parameter  characteristic  of  the  i-j  interaction  for  T 

c 

Parameter  characteristic  of  the  i-j  interaction  for  V 
Chemical  potential 
Acentric  factor 


Subscripts 

C  Critical  condition 

E  Excess  function 

r  Reduced  condition 

i ,  j , k, etc. Component  identification 
1,2  Component  identification 

ij  Interaction  between  i  and  j  species 
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Superscripts 

'  , '  ' , 1  '  ' , . Denotes  di fferent  phases 
~  Denotes  property  in  solution 

L  Liquid 

S  Saturated 

Vapor 


V 
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TABLE  1  :  Effect  of  Temperature  and  Pressure  on  the  Refractive  Index 
of  Carbonyl  Sulfide 


Pressure 

Refractive  Index* 

PSIA 

20°F 

60°F 

1 00°F 

1  40°F 

1 80°F 

76. 5+ 

1 .0062 

76 . 5++ 

1.3883 

113.0 

1.3877 

142. 5+ 

1.0098 

142. 5++ 

1 . 3666 

213.5 

1.3884 

1.3672 

245. 5+ 

1.0174 

245. 5++ 

1.3412 

313.5 

1.3891 

1.3682 

1.3424 

374. 5+ 

1 .0280 

374. 5++ 

1.3151 

413.5 

1.3894 

1.3691 

1 .3438 

1.3192 

513.5 

1.3898 

1.3698 

1.3451 

1.3209 

568. 5+ 

1.0444 

568. 5++ 

1.2863 

613.5 

1.3904 

1.3707 

1 .3464 

1.3227 

1.2881 

713.5 

1.3910 

1.3715 

1.3476 

1.3246 

1.2922 

813.5 

1.3916 

1.3723 

1.3487 

1.3267 

1.2953 

913.5 

1.3922 

1.3732 

1.3498 

1.3284 

1.2984 

1013.5 

1.3928 

1.3739 

1.3507 

1.3299 

1.3015 

JL  O 

Relative  to  vacuum  at  6328  A 
+  Saturated  Vapor 
++Saturated  Liquid 
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TABLE  2  :  Effect  of  Temperature  and  Pressure  on  the  Refractive  Index 
Of  Propane 


Pressure 

Refractive  Index* 

PSIA 

20°F 

60°F 

1 00°F 

1 40°F 

180°F 

55, 8+ 

1.0047 

55 . 8++ 

1.3128 

109. 5+ 

1.0088 

++ 

109.5 

1  .2938 

112.3 

1.3133 

190. 5+ 

1.0155 

190. 5++ 

1.2722 

213.5 

1.3142 

1.2949 

1.2727 

306. 5+ 

1.0271 

*j**f 

306.5 

1.2465 

313.5 

1.3148 

1 .2959 

1.2744 

413.5 

1.3156 

1.2970 

1.2757 

1.2498 

475. 5+ 

1.0491 

■H" 

475.5 

1 .2101 

513.5 

1.3163 

1.2980 

1.2773 

1.2525 

1.2131 

613.5 

1.3170 

1.2989 

1.2787 

1.2550 

1 .2206 

713.5 

1.3177 

1.2996 

1.2801 

1.2572 

1.2264 

813.5 

1.3183 

1.3006 

1.2814 

1  .2593 

1.2308 

913.5 

1.3189 

1.3015 

1.2830 

1.2613 

1.2340 

1013.5 

1.3195 

1.3023 

1.2838 

1 .2631 

1.2373 

* 

O 

Relative 

to  vacuum 

at  6328  A 

Saturated  Vapor 
++Saturated  Liquid 
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TABLE  3  :  Refractive  Indices  of  Carbonyl  Sulfide  at  Saturated  Conditions 
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TABLE  4  :  Experimental  Equilibrium  Phase  Properties  for  the  Carbonyl  Sulfide-Propane  System 
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TABLE 


5  :  Chueh-Prausni tz  Parameters  for  the  Pure  COS  and 
and  its  Binary 


COS 

=  0.056 

n 

av 

n 

=  0.33 

=  0.152 

=  0.49 

C3 

aL 

12 

12 

=  -48.71°R 
=  -0.353ft3/lb-mol e 

bv 

fik 

Dl 

=  0.075 

=  0.098 

i  J 

=  0.0044 

L 

'emperature  (°F) 


£] oClk^QWcu.ft. ) 


20 

60 

100 

140 

180 


0.0363 

0.0470 

0.0460 

0.0335 

0.0151 


TABLE  6  :  Comparison  of  the  Chueh-Prausnitz  Predictions  and  Experimental  Equilibrium  Ratio, 
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TABLE  7  :  Prediction  of  the  Critical  Properties  of  the  Carbonyl 
Sulfide-Propane  Binary 


impositi  on 

Temperature 

Pressure 

Molar  Volume 

xcos 

PSIA 

f t^/1 b-rnol e 

0.0 

206.20 

618.7 

3.180 

0.0500 

203.60 

606.3 

3.116 

0.1000 

201.28 

594.8 

3.052 

0.1500 

199.24 

584.4 

2.990 

0.2000 

197.51 

575.4 

2.928 

0.2500 

196.11 

568.0 

2.869 

0.3000 

195.06 

562.7 

2.810 

0.3500 

194.39 

559.6 

2.753 

0.4000 

194.12 

559.4 

.  2.698 

0.4500 

194.29 

562.3 

2.645 

0.5000 

194.91 

568.9 

2.594 

0.5500 

196.03 

579.7 

2.546 

0.6000 

197.69 

595.2 

2.499 

0.6500 

199.90 

616.1 

2.456 

0.7000 

202.73 

642.9 

2.415 

0.7500 

206.20 

676.2 

2.378 

0.8000 

210.36 

716.7 

2.344 

0.8500 

215.27 

764.8 

2.314 

0.9000 

220.97  ' 

'821.2 

2.288 

0.9500 

227.53 

886.2 

2.267 

1.0000 

235.00 

960.0 

2.250 

- 

* 

.  • 

TABLE  8  :Thermodynamic  Consistency  Test  of  the  Data* 
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FIG.  20  GAS  CHROMATOGRAPH  CALIBRATION  FOR  CARBONYL  SULFIDE-PROPANE 
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